changed so much by the overlap with the objects of such very small mass scales. 
Hence, as for the very large mass scale M2, the mass function of M2 itself might not 
be much affected by the effect of the overlap. The effect of the overlap on medium 
mass scales might give a serious problem for estimating the mass function. 

REFERENCES 

Bardeen, J. JVL, Bond, J. R., Kaiser, N. & Szalay, A. S., 1986, ApJ, 304, 15 

Bond, J. R., Cole, S., Efstathiou, G. & Kaiser, N., 1991, ApJ, 379, 440 

Jedamzik, K., 1994, preprint, astro-ph/9408080. 

Peacock, J. A. & Heavens, A. F., 1985, MNRAS, 217, 805 

Peacock, J. A. & Heavens, A. F., 1990, MNRAS, 243, 133 

Peebles, P.J.E., 1993, The Principles of Physical Cosmology 
(Princeton: Princeton Univ. Press) 

Press, W. H. & Schechter, P., 1974, ApJ, 187, 425 



23 



We assume that there is an isolated collapsed object of mass scale M2 overdense 
by just 8 C with a maximum peak density at r = 0. If an isolated object of mass 
scale Mi(< M 2 ) locates at R2 — R\ ^ r ^ R2 + i?i, the overlap between the objects 
of mass M2 and Mi occur. 

Now we consider the condition that the object of mass scale Mi, whose scale 
is between mass scale and M2, has the density of S c at R2 — R\ ^ r ^ R2 + R\ , 
that is, the object of mass scale Mi overlap with the object of mass scale M2. So 
we consider the following conditional probability of finding a region of mass scale 
between and M2 overdense by just S c at the distance R2 — R\ ^ r ^ R2 + R\ 
from the center of the isolated object of mass scale M2. 



when the mass scale M2 is much smaller than the characteristic scale M*, both 
conditional probability P(i?2 — R\ ^ f ^ R2 + Ri, 0, M2\peak) and P(i?2 — Ri ^ 
r ^ i?2 + M2, M2\peak) can be nearly equal to ^ because both <r(Mi) and 
<r(M2) are much larger than unity and then 7(r) is nearly equals to in each case. 
So the conditional probability P(6 = <5 C ,0 < M < M2) can be neglect for very 
small mass scales. But if the mass scale M2 is not much smaller than M*, the 
probability of the overlap might not be neglected. 

However, the objects with much smaller mass scale Mi(^ M2) dominantly 
contribute to the overlap with the object of mass scale M2 as we can see from the 
above probability. So the mass of the object with mass scale M2(^> M*) is not 



P(6 = <5 C ,0 < M < M 2 ) = {P(R 2 -Ri&r& R 2 + 0, M 2 \peak)- 

P(R 2 -Ri & r&R 2 + Ri,M 2 ,M 2 \peak)} 
I {I - P(R 2 -R!^r^R 2 + Ri, M 2 , M 2 \peak)}, 



where 



P(R 2 - Ri & r & R 2 + Ri , Mi , M 2 \peak) 




e 2 dy 
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s(r) = (1 + Ai(r) - 3cA 2 (r))a„(r) 



As in §4, we get the spatial averaged probability P(Mi, M2\peak) with the 
constant value defined as 

7 = , 1 (-)(! - s). 

V 1 -(fr^« 2 ( 1 + A 3) ^ 

s = (1 + Ai - 3cA 2 )a (^-13) 

Also in this case as in §4, we can approximate 7 as follows with good accu- 
racy(within 10% error around M\ ~ Af 2 ) because a ~ 0.5 and also we find that 
A 3 < 1: 

7~— (A - 14) 
We find that s is a little smaller than a. 

Hereafter, we can use the same procedure as in §4, in which the difference 
between the exact mass function derived by this formula and the PS mass function 
are deduced . Then we conclude that the deviation from the PS mass function is 
similar in the results shown in 54 while the factor s is a little different from a: 



sn ps (M M*) (A - 15) 

^s(l - s)(2 - s)(M- )l+f «5 c exp{i(f )(*+? )6l(l - sf}n ps (M » M* 



APPENDIX B 

In this appendix, we consider the effect of the overlap between the collapsed 
objects. 
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P(r,M 1 ,M 2 \peak) = 
1 

/Ml-^(1 + A 3 (r))) 



\2 



where 



, lc 1 2 l- £ g r (l + A 3 (r)) 111 J 



1 + Ai(r) = 1 + ^(n + 3)(n + 7)(1 - 

4 a„, ' 



1 , _ _ , a ra+ 2 • 



A2 ^V80 (n + 3)(n + 7)(n + 5)2(1 a. 
1 + A 3 (r) = 1 + ^(n + 3)(n + 7)(1 " ra+2 ^ 2 



4 a 



n 



v<5 + V7 + "8 = -3cz/ c 
Here a n is defined as follows: 



(7 



(0) J 1 1 &r 



where VKi = W(k,R\) and W2 = W(k,R-2). When the isolated object has the 
maximum peak density, v%,vj and v% must be negative value and so c must be 
positive. And we also find that its value is around 1 for any n. 

Furthermore, the conditional probability is rewritten by 



P(r,M 1 ,M 2 \peak) = —= I e~**dy, (A — 12) 

V27T J 7 ( r ) 

where 

v lc - e 0r (l + Ai(r) - 3cA 2 )z/ 2c «5 C 1 - a n (r)(l + Ai(r) - 3cA 2 ) 
7l-eg r (l + A 3 (r)) 

_ #c 1 ~ s(r) 

^ ^l-ejW(r)(l + A 3 (r)) 



^l-elal(r)(l + A 3 (r)) 
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multivariate Gaussians as follows: 



P(r,M 1 ,M 2 \peak) = -±— [°° P{\)dv u (A - 7) 

nv ) j vic 



where 

V = (5(r),5,$,^,5j-) 



V = (6,6-, 6^,6'^] 
From eqs.(A-l) and (A-7), we obtain 

P(r,M 1 ,M 2 \peak) = P{vi > v lc \v 2 , . . . , v n ) 



1 f°° — 1 O 

exp( — y±)dv u (A - 8) 



V^ 7r Q2 J Vlc 2 $2 



9£2i-£Q2 

^1 - £0r{ .92 J^ 2 + 7 _ 9 2 ^ 6 + " 7 + ^ 

1 ~~ 5 £ 02 ^ 1 — 5 £ 02 

2 
2 

2 



1 9£2i-£Q2 q 

5e 0r 1 <J £ 0r£02 — £2r / x-,2 



1 9 £2r 18g2r£Q2 



'Or q F 2 

1 "^02 

- 1 - 5 



where v l and ^ are defined as normalized variable, 



Vi 

v, = 



1 5 

2\2 



^lc = — 



W> ! (4-9) 



r 



And e is correlation coefficient, 

e r = e(r) 



£02 



(70Cr r 

3 (70(72 
3 <J 2 (7 r 



, 2r = _^ji_ ( ,4- 10) 



From the equations (A-6)~(A-f 0), we can rewrite the conditional probability 
as follows: 
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M 22 = al 

2 

M 33 = M 44 \/-,-, .,' 

M 66 = M 77 = M 88 = -± 

5 



i 

15 



M99 = .W|„.|„ = ,\/||.n = M 67 = M 78 = M 86 - n - 

2 

M 26 = M 27 = M 28 = CTl 



M 2 i = (rl r 



3 



r 2 



W|„ W|; \/|. .J'. (A -5) 



where 



4ttV 

^ f27r r ./o 

c 00 



i = ^l W*(k,R 2 )\6 k \W m +Qdk 



a r = 7^ W 2 (k,R 1 )\S k \ 2 k 2 dk, (A-6) 
\ ZlT ) Jo 

and W(k, R) is a window function and we find that a 2 r = a^a(r) which is defined 
in §4. Throughout this paper, we use the sharp-A; filter function, therefore, W(k } R) 
is given by 

f 1 (k < MR)) 
W(k,R)= \ y ~ y " 
1 (k > k c (R)) 

By calculating these components, we obtain the multivariate Gaussian distribution 
function P(\). 



(b) Conditional probability of maximum peak 



Now, the probability which we want is the conditional probability 
P(r, Mi, M 2 \peak) (see §5). From Bayes' theorem, this probability is related to 
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mass function approximately. 



(a) Multivariate Gaussian distribution function 
When m-variables are associated with a Gaussian field, the probabilities given 

by 



(-#) 

P(V)dV = , 6 2 =dV, (A -I) 

V ' ^(27r)™det(M) V 



where 



q = v t m _1 v (A -2) 

V = (yi,...,y m ) (A-3) 
Mij = ((yi — (yi))(yj — (yj))} (A - 4) 



We use the brackets () as the ensemble mean of the universe, but in practice, 
assuming homogeneity and ergodicity in space , we take it as the spatial mean, so 
(yi) corresponds to the spatial mean of y t and Mij is covariance between y t and 
yj. Now, because we want the conditional probability of maximum peak, we must 
calculate the first and second derivative of the density perturbation 8. So we need 
to consider the 11 X 11 covariance matrix M, where V = (<5(r), 8, S' l} S'- l} 8'-j), z, j = 
1 , 2, 3, z 7^ j, and ' means derivative with respect to r. The non-zero components 
of this matrix follows (of course, M is symmetric matrix): 
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to estimate the mass function when we consider the effect of the spatial correlation 
in the cloud-in-cloud problem. By using this formalism, we showed how the effect 
of the spatial correlation alter the PS mass function. Furthermore, we showed 
the exact formula of finding the mass function in taking into account the exact 
probability that the isolated object has the maximum peak density. Although it 
is difficult to solve numerically mass function with this formula in practice, we 
can approximately show that the deviation of the modified mass function from 
the PS mass function is similar in the case that we consider only the effect of 
the spatial correlation. We conclude that while the deviation on smaller mass 
scales is within factor a, whose value depends on the power spectrum, it is very 
large on larger mass scales in some cases. So we warn us that we must be careful 
about the limitation of the PS formalism when we apply the PS mass function 
to some problems. Furthermore we found that we can neglect the probability of 
overlap between the collapsed objects on very small mass scales while it might 
not be neglected on other mass scales. This problem of the overlap especially on 
medium mass scales is a serious one to estimate the mass function in the statistical 
argument such as the PS formalism or our formalism. 
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APPENDIX A 

In §5, we mentioned the exact formula of mass function. Here, we express 
the conditional probability on which the exact formula is based. And we will also 
show how the mass function derived by the exact formula deviate from the PS 
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The derivation of the above probability is explained in Appendix A. Thus we 
get the conditional probability P(Mi, M2|peak) after averaging spatially the above 
probability. This is the exact procedure of getting the mass function. In practice, 
however, it is much complicated to solve numerically the mass function. But we can 
see how the mass function with this formalism deviate from the PS mass function 
approximately. The detailed method is explained in the Appendix A. The result 
is that the deviation is similar to that shown in §4(eq.l7), in which we neglect 
the condition that the object with mass M2 has maximum peak density. So we 
conclude that the deviation is within factor on smaller mass scales, however, it is 
very large on larger mass scales in some cases(see eq.(A-15)). 

Furthermore we investigate the probability that the collapsed objects would 
overlap. If the overlap would occur, we could not account the number of the objects 
because we can not determine whether the overlapped objects merge to one larger 
mass object or fragment into two small objects in our statistical argument. In this 
case it is necessary to analyze dynamically the effect of the overlap. As shown 
in the Appendix B, by using the exact formula of mass function, we find that 
the probability of the overlap can be neglected on very small mass scales while it 
might not be neglected on other mass scales. The effect of the overlap especially 
on medium mass scales gives a serious problem to estimate the mass function in 
the statistical argument. 

6. CONCLUSIONS AND DISCUSSIONS 

We have reanalyzed the PS formalism. When we neglect the spatial correlation 
of the density fluctuations in the cloud-in-cloud problem, we can get just the PS 
mass function including the factor of two in the sharp k-sp&ce filter by using the 
Jedamzik formalism, which is different approach from the previous one used by 
Peacock & Heavens(I990) and Bond et al.(1991). However we believe that it is 
very important to take into account the effect of the spatial correlation in order to 
estimate correctly the mass function in real. The Jedamzik formalism is very useful 
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So we need to consider the conditional probability of finding a region of mass scale 
Mi overdense by 8 C or more at distance r, provided the object with mass scale is 
included in an isolated object of mass scale M2 with the maximum peak density of 
&M 2 = <*> c at r = 0. Then we must use the following constrained probability, 

P(r, Mi, M 2 (peak) = P(6 Ml (r) > 6 C \6 M2 = <5 c ,peak) 

f f 00 ,-\ 



',2 Jv lc Z V2 



1 9g2r£Q2 q 

^ r r 1 5e 0r I , 6 £ 0r£02 ~ £2r , . . 

h = [Vi - £0r{ .92 J^ 2 + 7 i 9 2 ^ 6 + ^ + ^8jJ , 

1 ~~ 5 £ 02 ^ 1 — 5 £ 02 

2 

_|_ 9 ^2r 18g2r£Q2 

/2 - J- - £Qr 



q F 2 5 

1 "^02 

1 5 

where v l and z/i c are defined as normalized variables: 



Of) 1 ' 



I 
Oc 

Vic = — , 
(T r 

V = (<5(r),<5,(5-,^,^). 
And e is correlation coefficient, 



£0r 



£02 



£ 2r 



(70Cr r 

3 (70(72 
3 (72 (7 r 



where cr means variance, 
4ttV 



a 2 





00 



W 2 (k,R 2 )\S k \ 2 k^ 2m+2 Uk 

./n ftf 



^ = 7^3 / W 2 (k,R 1 )\6 k \ 2 k 2 dk. 
{ Z7r ) Jo 
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In this case, 6 c /cr(M) ^> 1. We can use the next approximate expression, 



1 1 - 

2 dx ~ e 2 ^ (/3 > r 



2tt J/? V2tt P 



So, the ratio i 



is 



J r f , ,2df d r /\ ,2df 



dM g dM dx g dx 



-y/2w(l ~ a){l - (1 - a) 2 }xexp{-x 2 (l - a) 2 }, (16) 



where x = S c /a = (M/M*)5 + %6 C 



From the above equations (15) and (16), the modified mass function is given 
as follows: 

an ps (M < M*) (17) 

&a{\ - a)(2 - «)(f )^^ c exp{i(f )d+f )«5 C 2 (1 - af}n ps (M » M*) 

As we can see, the deviation from PS mass function n ps is within factor a 
on smaller mass scales(Af <C M*)- However on larger mass scale ranges of mass 
function(Af ^> M*), we find large deviation from the PS mass function in some 
cases. Of course, it must be noted that the characteristic mass M* change as time 
goes while M* is initially small. 



5. EXACT FORMULA OF MASS FUNCTION 

In the previous section, we showed how the effect of the spatial correlation alter 
the PS mass function. Strictly speaking, the conditional probability P(M\, M2) 
shown in §4 is not sufficient to get the exact mass function. The isolated object 
with a large mass M2 must have the maximum peak density with 8m 2 = <*> c at r = 0. 
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Using the Jedamzik formalism(eq.(5)), 



f(>*c,M)_i r~ d P(M,M>) 

Jo { ' g(M) { ' 



g(M) p Jo 

Here, 

P(M, M') J 1 (M <M') 
g{M) ~ \ (M >M') 

Then we can get the following equation, 



d -{ f -}\ = -Mn(M) (14), 



'dMg JI p 



because d{f}/dM = -S(M - M'). 

Now we estimate the mass function in two extreme cases. First of them, we 
consider the small mass scale ranges of mass function (M <C M*). In this case, we 
can use the relation such as 6 c /cr(M) <C 1. Then 

f = ~lx= / & 2 dy 



2tt J -4^ 2 ^o(M) 

ct(M) v v ' 

1 r °° - y ~, 1 1 S ° n ^ 
e 2 dy ~ -= (1 — a) 



In order to see the difference between the modified mass function and the PS mass 
function n pS} we see the ratio of d(f / g) / dM with 2df/dM. 



dM X g ]l dM lf ~ 1/2 2df lf ~ 1/2 



a/2 

2{(l-a)/ + a/2}^= 1 /2- a ( 15 ) 



Next we consider the larger mass scale ranges of the mass function (M ^> M* 
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fined as 

«R i-R 



f It fit 

I'i \/|. / I'ir. M\. M. )\-r-<lr / !-/-'<//• 

Jo Jo 

e'^dy, (10) 



1 r°° .,2 



V27T 

where i? ~ i? 2 = (3M 2 /4tt / o) 1 / 3 and 

= , 1 -(1 - a) = . 1 = (%! - a). 

i - ei& ^ Twfr^ ai 

The above equations give the definitions of the constant values, f3 and a. Then 
after when we calculate this conditional probability in any case, we can solve the 
mass function by substituting this probability into the integral equation (5). If we 
neglect the spatial correlation, then a = a(0) = 1 and so P(Mi, Af 2 ) = 1/2. Hence 
again we can get the PS mass function. 

In practice, it may be a little complicated to solve numerically the mass function 
with the integral equation. However we can get approximately the mass function 
with good accuracy as shown below: we find that a is between and 1 and its 
value depends on the power spectrum, but a is approximately around 0.5 (for 
—2 ^ n ^ 2). We consider these cases for that a 2 is much less than unity. In these 
cases, we can approximately get 

(1-a). (11) 

o"l 

The error at M\ ~ Af 2 is at most around 10%. Then we get the following condi- 
tional probability with good accuracy: 



P(M, M') \ " ^ = 7T« If e " V ^ ( M ^ M ') 
= (M > M') (12) 



where f3 = ^y(l - a). 
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correlation, while in this probability the condition that the 8m 2 has the maximum 
peak density is neglected(see §5). 

We find that P(r, Mi, M2) is equivalent to the constrained probability defined 
as follows (Bardeen et al. 1986): 

P(r,M u M 2 ) = P(S Ml (r) > 6 C \6 M2 = 8 C ) 

exp[ — — 7— — - — \dv\ (6) 



sir) = eoa(r) = -^—a(r), (8) 



^(l-e 2 (r))J Vlc 2(1 - £ 2 (r) ) 

where v\ , V2 have a unity variance define as 

v\ = — , v 2 = — , cri = <t(Mi), ct 2 = cr(M 2 ), = — , ^2c = — (7) 

G\ (72 (71 (72 

and e(r) is defined by 

i 

(71 (72 

which corresponds to the two-point correlation function. Here cr/j is a cross- 
correlation of the density fluctuations with the mass scale M\ and the mass scale 
M2 at the same point. When we use a sharp &-space filter, the cross-correlation 
is the same as the variance of the density fluctuation filtered with the larger mass 
scale. Therefore, is equal to the variance 02 ■ Then, in this case, £0 = 02 /(T\. 

The conditional probability is rewritten by 

P(r,Mi,M 2 ) = -= / ( ' <///. (9) 

V27T Jp{r) 



where 



0/ s v lc -e(r)v 2 c $c 1 - a(r) 



Vl-e 2 (0 CT i y/l- e 2 a 2 (r) 
v\ - e a(r)v2c 



1 — £o« 2 (r) 



Now we consider the spatial averaged conditional probability P(Mi, M2) de- 
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the sharp k-space filter, the factor of 2 is correct even when we use the Jedamzik 
formalism, which is different approach from the previous works (in Jedamzik's 
paper(f994), he used the wrong conditional probability in the sharp k-space filter. 
So he derived the wrong result that the factor of 2 is incorrect even using the sharp 
k-space filter). This Jedamzik formula is very useful to analyze the mass function 
with the effect of the spatial correlation functions in the cloud- in-cloud problem. If 
we can get the conditional probability P(M, M') in taking into account the spatial 
correlation, we can in principle estimate the mass function by solving the integral 
equation (eq.(5)). 



4. SPATIAL CORRELATION 

In the previous works, they considered the probability for the density fluctua- 
tions only at one point. In reality, however, the object of mass scale M will have 
a finite size and will include the other objects in this finite region. So we must 
consider the two-point correlation of the density fluctuations in order to analyze 
whether the objects of smaller mass scale Mi are included in the finite region of 
the object with a larger mass scale Mi- Then we believe that it is necessary to 
consider the spatial correlation of the density fluctuations in order to solve the 
cloud-in-cloud problem correctly. 

Now we consider the following conditional probability for simplicity in order 
to analyze only the effects of the spatial correlation on the PS formalism at first (in 
the next section, we will show the exact conditional probability for estimating the 
correct mass function while the form of the probability is complicated). In esti- 
mating the conditional probability P(Mi, M2), at first, we consider the conditional 
probability, P(r, M\ , M2), of finding a region of mass scale M\ overdense by 8 C or 
more at distance r from the center (r = 0) of the isolated object of mass scale M2, 
provided the object with mass M\ is included in a object of mass scale M2(> Mi) 
with a finite size and 8m 2 = <*> c at r = 0. As we will show later, it is enough good 
to consider this conditional probability in order to analyze the effects of the spatial 
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Here n(M) is the number density(«.e. mass function) of the isolated objects of 
mass scale M. If we can calculate P(Mi, M2) in any case, we can get a mass 
function n(M) without worrying about the cloud-in-cloud problem in principle by 
using this formula(we call it the Jedamzik formalism). 

By definition, the conditional probability must satisfy P(M } M') = for M > 
M'. As if P(M, M') would have the property such as P(M, M') 1 lor M 
then 

<//' M 

1 J 1 —n(M), 



'dM 1 p 

where we use the relation such as dP(M } M')/dM = —6(M — M'). Then we get 

n{M)=*-\*L\. 
v ; M d.\l 

This is the original PS formula without the factor of 2. In general, however, this 
case could not appear in real. 

In the sharp &-space filter, we find 

P(M,M')= f%-^l exp{-i^^} 
Js c V^^sub ^ a gnh 

_ 1 

" 2' 

where 

CT 2 _ f ^ r (l + f)_ f ^ r (l + f) 

swh M. M. 

Here we use the fact that the region with mass M' is the isolated one, then Sj^i = 8 C . 

Then dl'i \l.\l' \ d.\l = -\6{M - M'). Hence, from eq.(5), 

1 J 1 —n(M) 



'dM 1 2 p 
and so 

P 1 df 



n(M) = 2- 



M dM 

This is just the PS mass function including the factor of 2. Then we proved that in 
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Press and Schechter simply multiply the probability by a factor 2 with no clear 
physical reason. Thus, 

n(M) = 2^-\4L\ = rips- 

The factor of 2 has long been noted as weak point in the PS formalism. The 
problem of correctly counting the underdense regions which are embedded with 
overdense regions is called "cloud-in-cloud" problem. Peacock & Heavens(1990) 
and Bond et al.(1991) proposed a solution to the cloud-in-cloud problem in taking 
into account the probability, P up} that subsequent filtering with larger scales might 
at some points result in having 8 > 8 C when at smaller filter 8 < 8 C . Surprisingly, 
in the sharp k-space filter, it is found that the factor of 2 in the PS formalism is 
correct. However, the factor of 2 is incorrect when using the other filters. 

3. JEDAMZIK FORMALISM 

Jedamzik(1994) proposed the another approach to the cloud-in-cloud problem. 
We will explain this approach below: our goal is to compute the number density of 
isolated region. Isolated regions are regions which have collapsed but not included 
in yet larger regions overdense by 8 C or more. In particular, a region of mass scale 
M will be only counted as an eventually vilialized object of mass M, if for any 
larger mass scale the average overdense of the region is below the critical density. 
Then the density fluctuation of the isolated regions is just the critical density 8 C 
because the regions with 8 > 8 C would be eventually counted as an object of the 
larger scale. 

Now we define the conditional probability P(Mi, M 2 ) of finding a region of 
mass scale M\ overdense by 8 C or more, provided it is included in an isolated over 
dense region of mass scale M 2 (> Mi). Then the volume ratio /(> <5 c ,Mi), of 
finding a region of mass scale Mi overdense by 8 C or more is given by the following 
integral form: 

f(>6 c ,M 1 ) = f°° dM 2 n(M 2 )^P(M u M 2 ) (5) 
Jo P 
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The probability /(> 8 C} M) to find a region of mass scale M to be overdense by 
more than the critical amount 8 C is given by an integral over the tail of a Gaussian 
distribution function, 

,00 1 S 2 

f(>6 c ,M)= exp{ — — }d6, (2) 

J{ ' } J Sc ^a 2 (M) VX 2a 2 (M) i ' 1 } 

because linearly estimated 8 is the random Gaussian fields. Here the variance 
a 2 (M) is estimated by sum up a variance of each Fourier components in the sharp 
A;- space filter as follows: 



v M M ) 

a2 ( M ) = 7^ \8 k \ 2 47rk 2 dk, (3) 

^7T ) JO 



k c (M) = ,l(^-_fl\ 

47T/J 



When we assume a simple power law, P(k) oc k n } the variance a (M) is given by 

o M + 3 ) 

AM) = (-)-r- (4) 

where M* is the characteristic mass scale such that cr 2 (M*) = 1. 

In the PS formalism, /(> 8 C , M) is considered to be proportional to the prob- 
ability that a given point has ever been processed through a collapsed objects of 
mass scales > M. If 8m > 8 C for a given mass M at a point, then it will have 
8 = 8 C when filtered on some larger mass scale and so it will be counted as collapsed 
objects of the larger mass scale. So in this argument, it is assumed that the only 
objects which exist at a given epoch are those which have just collapsed^ = 8 C ). 
Thus the mass function n(M) (defined such that n(M)dM is the comoving number 
density of the collapsed objects in the range dM) is given by 

— n( M = - — = , — exp 8 r a . 

In this argument, it only associates half of the mass of the initial density fluctua- 
tions with eventually to be collapsed. The problem is that half of the mass, which 
is initially present in underdense regions, remains unaccounted for. 
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fluctuations, is given by(Bardeen et al. 1986) 



P(|£kUkMIWk 



2\8 k 



exp{ 



}d\8 k 



#k 



(1) 



P(k) 



P(k) 



2tt ' 



where </>k is the random phase of 8^ and P(k) is the power spectrum. 



In general, we assume that the density fluctuations associated with the objects 
of mass scale M is given by the filtered density fluctuations defined as follows: 



where W(k } R) is an window function and p is the mean density. 

It is called "sharp k-space filter" when we use the window function defined as 



In this paper, we consider only this sharp k-space filter for simplicity. If we 
choose other filters, the results might change quantitatively. 

The density fluctuation 8m grows up as time goes and the objects with mass 
scale M recollapse and vilialize when the linearly estimated 8m grows to the critical 
density 8 C . It is well known that the critical density 8 C ~ 1.69 for the spherical 
collapse in the flat universe(Peebles 1993). Throughout this paper, we consider 
only the spherical symmetric case. Also in other cases, the critical density has the 
similar value. It is considered that the objects with mass scale M vilialized at this 
epoch. 
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that the factor of two in the PS formalism is correct in the sharp k-space filter if 
the effects of the spatial correlation is neglected. 

In this paper, we present how the effect of the spatial correlation in the cloud- 
in-cloud problem alter the PS formalism. And we show the new formula of mass 
function in taking into account the spatial correlation. In fact, it is very compli- 
cated to derive numerically the mass function with the new formula. And the form 
of the numerically derived mass function may not be simple enough for application 
to the various problems. However, we believe that it is very important task to 
analyze the limitation of the PS mass function. So we show approximately how 
the PS mass function deviate from the modified one in some special cases. 

As for the derivation of mass function, we use the Jedamzik formalism with 
the correct conditional probability including the effect of the spatial correlation. 

In §2, the PS formalism is briefly reviewed. And in §3, the Jedamzik formalism 
is shown. We present the modified formula with the spatial correlation in §4 and 
present how this effect of the spatial correlation alter the PS mass function. In §5 
we also show the new formula of mass function using the constrained probability 
that the isolated object has maximum peak density, whose condition is necessary 
in order to correctly get the mass function in taking into account the spatial corre- 
lation. In §6, we devote to the conclusions and discussions. In the Appendix A, we 
show how the mass function derived by the new formula deviate from the PS mass 
function approximately and also suggest the troublesome effect of the overlap of 
the collapsed objects in the Appendix B. 

2. THE PS FORMALISM 

In this section, we briefly review the PS formalism. In the standard cosmo- 
logical models, we assume that the initial small density fluctuations are the seeds 
of cosmic structures such as galaxies and they obey the random Gaussian dis- 
tributions. Then the probability function for the Fourier components of density 
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underdense regions which are embedded within overdense regions (it is called the 
"cloud-in-cloud" problem). 

Peacock and Heavens (1990) and Bond et al. (1991) attempted to solve the 
cloud-in-cloud problem by using the peak formalism (Peacock and Heavens 1985; 
Bardeen et al. 1986). They considered the behavior of values of density peaks under 
varying the length of filters for smoothing the density field, that is, the probability 
that the underdense regions upcross to overdense regions in different length filters. 
In the sharp k-space filter, because of random phase, the behavior of the density 
field when summed up a wave with larger wavenumber k } is expressed by random 
walk. Only in this case, surprisingly, they obtained the same result as the PS 
formalism including the factor of two. In other filters instead of the sharp k-space 
filter, it is found that the factor 2 introduced in the PS formalism is not correct. 
As for as the cloud-in-cloud problem, these analyses are, however, insufficient : the 
spatial correlation functions of density fluctuations must be considered in order 
to solve correctly the cloud-in-cloud problem. In the previous works, however, 
they considered only the probability for the density fluctuations at one point in 
the space and neglected the spatial correlation of the density fluctuations in the 
cloud-in-cloud problem. In fact, it is found that the effect of the spatial correlation 
alter the PS formalism even in using the sharp k-space filter as we will show later. 

Jedamzik (1994) approached the cloud-in-cloud problem with the formula using 
the integral equation of the mass function, which is different approach from the 
previous works. We believe it is useful approach to the cloud-in-cloud problem 
in taking into account the spatial correlation. But in his paper, he neglected the 
spatial correlation. Furthermore he used the wrong conditional probability about 
the condition of the isolated objects. So he misleadingly concluded that the factor 
of two in the PS formalism is incorrect even in the sharp k-space filter, which is 
inconsistent with Peacock & Heavens(1990) and Bond et a/.(1991). As we will 
show later, however, we find that the factor of two can be derived in the sharp 
A;- space filter even when we use the Jedamzik formalism in neglecting the spatial 
correlation, which is consistent with the previous works. Then we will conclude 
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1. INTRODUCTION 

When we would like to understand the universe, the problem of the structure 
formation is one of the most important problems. In the standard cosmological 
model, we believe that the seeds of the structures are small density fluctuations in 
the early universe, and after the epoch of matter-radiation decoupling at redshift 
z ~ 1000, they have evolved and hierarchically clustered to the bound objects with 
extended magnitude of order of mass such as galaxies, clusters of galaxies and 
Lyman-a clouds, etc., via gravitational instability. 

Therefore if the mass function for these bound objects would be known, we 
can in principle obtain the informations of the nature of the initial density fluctua- 
tions. It is, however, rather difficult to obtain completely analytic solutions for the 
problem of the structure formations since it is too complex to treat the non-linear 
phases of collapse to bound objects. 

One of approaches to this problem, therefore, is N-body simulation. Since 
explicitly calculating gravitational force, it can follow the non-linear phases. How- 
ever, by present computers, N-body simulations can work only in a limited dynamic 
range of mass. 

Another analytic approach has been developed by Press and Schechter (1974; 
hereafter PS). They assumed that the initial density fields are expressed by the 
random Gaussian, and the bound objects have hierarchically formed from the initial 
density fluctuations . Then they proposed the formalism for counting the number 
density of objects with mass M (i.e., the mass function). The PS mass function 
has been applied to various problems such as the number densities of galaxies 
and clusters of galaxies because of its simple form. However, the problem of their 
argument is that half of the mass, which is initially present in underdense regions, 
remains unaccounted for. With no clear reasons they simply multiplied the number 
density by a factor of two to take into account underdense regions since all of the 
mass must be considered. So it has been recognized that this factor of two is a 
crucial weak point of the PS analysis. To solve this problem, we must consider 
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ABSTRACT 

The Press-Schechter(PS) formalism for the mass function of the collapsed ob- 
jects are reanalyzed. The factor of two in the Press-Schechter formalism is con- 
vinced to be correct in the sharp k-sp&ce filter even when we use the another 
approach proposed by Jedamzik(1994) in the cloud- in- cloud problem, which is 
different from the previous approach by Peacock & Heavens(1990) and Bond et 
al.(1991). The spatial correlation of the density fluctuations, however, had been 
neglected in the cloud-in-cloud problem. The effects of this spatial correlation is 
analyzed by using the Jedamzik formalism and it is found that this effect alter the 
PS mass function especially on larger mass scales. Furthermore the exact formula 
of deriving mass function is shown. We also find that the probability of the overlap 
of the collapsed objects can be neglected on very small mass scales while it might 
not be neglected on other mass scales. 

Subject headings: cosmology:theory-mass function- large scale structures-galaxy for- 
mations 



